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Abstract 



We discuss the design and performance of a threshold Cerenkov counter 
for identification of charged hadrons. The radiator is pressurized gas, which 
is contained in thin-wahed cyhndrical modules. A mirror system of novel 
design transports Cerenkov photons to photomultiplier tubes. This system 
is compact, contains relatively little material, and has a large fraction of 
active volume. A prototype of a module designed for the proposed CLEO III 
detector has been studied using cosmic rays. Results from these studies show 
good agreement with a detailed Monte Carlo simulation of the module and 
indicate that it should achieve separation of pions and kaons at the 2.5cr-3.0cr 
level in the momentum range 0.8-2.8 GeV/c. We predict performance for 
specific physics analyses using a GEANT-based simulation package. 



I. INTRODUCTION 

The identification of stable charged hadrons at momenta above 1 GeV/c is difficult to 
accomplish, particularly in a colliding beam environment where spatial limitations constrain 
detector design. This problem has received increased attention in the past few years in the 
context of flavor-dependent processes (e.g. CP violation) in decays of B mesons, where pions 
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and kaons must be identified with high efficiency and accuracy. To study such processes in 
e'^e~ colhsions at the T(4S) resonance will require devices which can provide identification 
for particles with momenta up to ~ 4 GeV/c. 

Presented here is a design for a particle identification (PID) device based on the Cerenkov 
effect, working mainly in the threshold mode with a radiator of pressurized gas. The detec- 
tor concept is generally applicable in situations with moderately high rates and low particle 
multiplicities. It has been applied here to the proposed CLEO III detector [Q], an upgrade to 
the CLEO II detector at the Cornell Electron Storage Ring, and some of the detector param- 
eters reflect our accommodation of constraints imposed by the existing system. Speciflcally, 
the radial gap available for a PID device was flxed at 20 cm, between the main tracking 
chamber (at a radius of 82 cm) and the calorimeter. In addition, it was prescribed that the 
PID device should comprise no more than 0.12 radiation lengths of material, so as to have 
a minimal impact on the performance of the Csl calorimeter. We flrst discuss the overall 
approach used to arrive at a design, including development of components and tests needed 
to establish validity. A description of a working prototype and tests with cosmic rays follows. 
Finally, results from these tests are used to tune a Monte Carlo simulation, which is then 
used to evaluate the performance of a full-scale system in the proposed CLEO III detector 

for physics analyses in symmetric e~^e~ collisions in the T region {^/s ~10 GeV). 

II. DETECTOR DESIGN 

A threshold Cerenkov counter has several advantages over other PID devices. Its readout 
can be simple and the number of channels few. The response time is intrinsically fast, and 
timing and pulse height calibration requirements are minimal. These features make the 
offiine software requirements rather simple. Furthermore, such a device may be included in 
the trigger. 

Design considerations are dominated by the need to maximize photon detection efficien- 
cies, because the absolute number of photons produced by Cerenkov radiation is relatively 
small. It is also desirable to restrict the material thickness of the detector and to minimize 
inactive volume within the central region. Other considerations include ability to operate 
in high magnetic flelds, robustness, cost, and personnel requirements. 

As is described below, the index of refraction required for the desired threshold limits 
the choice of radiator to pressurized gases or silica Aerogel [^. We pursued the possibility 
of pressurized gas because it has better transparency, particularly at shorter wavelengths, 
where the Cerenkov spectrum is concentrated. Another advantage of gas is that there is 
some flexibility to vary the threshold through adjustments of the gas pressure. The high 
transparency enables the design of a long cylindrical detector module, in which radiated light 
is transported up to 2 meters before being detected. The modules may thus be arranged to 
span the central solid angle of a collider detector, and are sufficiently large that the number 
of readout channels required is small, of order 100. This system can be made conservatively, 
relying only on readily available components and proven technologies. More speciflc aspects 
of the design are discussed below. 
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A. Radiator 



In many detector systems, and specifically in this case, hadrons with momenta below 
0.8 GeV/c are readily identified through measurements of specific ionization (dE/dx) in 
tracking devices, so that the desired range of operation for an additional PID device is in 
the region above 0.8 GeV/c. To provide good identification at 0.8 GeV/c we select an index 
n = 1.02, which places the pion threshold at 0.70 GeV/c. The kaon threshold for this index 
is 2.45 GeV/c, which is near the kinematic limit for products of B decay in the rest frame 
of T(4S) events. 

For Cerenkov photons to be detected, the radiator must transmit them, and the photode- 
tector must be sensitive to them. For a system sensitive to photons with energies between 
Emin and Emax (in eV), a particle traveling with speed v = Pc a distance L (in cm) through 
a medium with refractive index n produces 

^ 370 sin^ ecL{E^ax - E^in) (1) 

detectable photons, where 6c = cos^^(l//5n) is the Cerenkov angle. The most effective 
device for detection of individual photons is still the conventional photomultiplier tube, 
which maintains a high quantum efficiency for wavelengths below 600 nm and may detect 
photons with wavelengths below 200 nm if the window is made of a medium transparent 
to them, such as quartz. Sensitivity to shorter wavelengths is desirable, as the Cerenkov 
spectrum is more concentrated in that region. However, because many common impurities 
absorb strongly at wavelengths below 200 nm, the practical range for detection is 200- 
600 nm (2.1-6.2 eV). In this range, in a material with index n = 1.02, 1 GeV/c pions radiate 
30 Cerenkov photons/cm, while 2.8 GeV/c pions radiate 60 photons/cm. If the detector 
system is able to count the number of photons, some separation of pions and kaons can be 
maintained even above the kaon threshold; only 11 photons/cm are emitted by 2.8 GeV/c 
kaons. Thus, pions and kaons at momenta in the range 0.8-2.8 GeV/c produce Cerenkov 
signals which should be readily distinguishable from each other. 

The selection of Cerenkov gases at our chosen index is limited, as many either liquefy 
under pressure at room temperature before reaching the target index or require undesirably 
high pressures. Sulfur Hexafluoride (SFg) satisfies our criteria. At an absolute pressure 
of 20 atmospheres, its room-temperature index of refraction is n = 1.02. Its inertness, 
non-fiammability, and previous use as radiator make it an attractive candidate for 
this application. Ethane (C2H6) is an alternative which reaches the same index at some- 
what higher pressure but has a longer radiation length. The higher pressure required, in 
combination with its fiammability, dissuaded us from pursuing ethane for this design. 

For pathlengths up to 1 atm-m (1 m pathlength at 1 atm), SFg is known to have excellent 
transparency to light with wavelengths above 190 nm Q]. However, because this application 
requires good transparency over a pathlength of 20 atm-m, we performed an independent 
measurement using a McPherson 218 0.3 m scanning monochromator illuminated by a deu- 
terium lamp. The monochromator was made to scan through wavelengths from 200 nm 
to 350 nm. The resulting beam traversed 1 meter of gas contained in a pressure vessel 
with sapphire end windows and was detected in an EMI-Gencom RFI/B-214 solar-blind 
photomultiplier tube. The anode current of the photomultiplier, which is proportional to 
the intensity of the transmitted light, was converted to a DC voltage by a Keithley 440 
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picoammeter, digitized by an AR-B3001 ADC plug-in card (Acrosser Technologies, Inc.), 
and read out at 0.25 nm intervals. The wavelength scanning and data acquisition of this 
system were controlled by a personal computer (IBM AT286 equivalent). 

Two grades of SFg were studied, instrument (99.99% pure) and ultra-pure (99.996%). 
As a baseline reference we used ultra-pure nitrogen (99.9999%), which has a transparency 
comparable to vacuum. For a pathlength of 1 atm-m, both grades of SFe were found to be as 
transparent as nitrogen above 200 nm, in agreement with [Q. Each gas scan was therefore 
normalized to its own 1 atm baseline rather than to nitrogen, reducing substantially the 
systematic error from variables such as changes in alignment between gas fills or drifting of 
electronic calibrations over time. 

The gas vessel was initially filled with SFg at the maximum pressure delivered by the 
cylinder (~20 atm) and allowed to stabilize [Q. Scans were then taken at different pressures 
by releasing gas successively until the ambient pressure was reached, at which point a scan 
was taken to act as the baseline for that set of measurements. The resulting transmission 
efficiencies as a function of wavelength are shown in Figure |I[ 20 atm-m of both instrument 
and ultrapure SFg were found to have transparencies comparable to 1 atm-m nitrogen, down 
to 320 nm. However, at shorter wavelengths, the instrument grade is more transparent than 
the ultra-pure, confirming the results in [Q. More than 90% of light down to 240 nm is 
transmitted by the instrument grade, dropping to 72% by 200 nm as Rayleigh scattering 
becomes more prevalent. Although it is not known to this group why the lower purity gives 
better performance, it has the fortuitous result of greatly reducing the cost of the gas needed 
for this detector. 

B. Pressure module 

The pressure vessel is designed with a cylindrical shape to insure the required mechanical 
strength with minimum wall thickness. Within each module, photons are transported via 
refiections off mirrored surfaces to the tube ends, where they enter a transparent window and 
are collected by photomultipliers. The vessel itself should be composed of a strong, low-Z 
material to minimize the number of radiation lengths of material penetrated by particles. 

In designing a device for CLEO III we envision modules with ~10 cm diameter and 2 m 
length. Two close packed layers, arranged in a barrel, surround the interaction point just 
outside the central tracking chamber (Figure ^). An economical option for the tube material 
is aluminum. A 0.1 cm wall thickness is within the limits prescribed for CLEO III and is 
adequate from the point of view of safety. The hoop stress for this design is about 1 x 10^ Pa, 
while Al tubes of grade 7075-T6 will maintain their strength to about 5x10* Pa, providing a 
five-fold safety factor. Other materials considered include carbon fiber and beryllium. These 
could reduce the material thickness but were not studied fully. 

C. Optics 

The relatively small diameter and long transport distance of the envisioned module place 
severe constraints on the optical design. To obtain good collection efficiency, photons should 
experience no more than 4-5 refiections before reaching the photomultiplier, as refiectivities 



4 



in the ultraviolet region can be expected to be ~ 80% at best. The simple option of 
placing a smooth cylindrical reflecting liner in the pressure vessel works remarkably well 
for tracks which are moderately angled toward the tube ends, at polar angles of 65° or less 
(|cos^| > 0.42). We used a mylar sheet, 0.076 mm thick and coated with a thin fllm of Al 
and a MgF2 overcoating |]^, which was rolled and placed inside the tube. Measurements on 
one sample indicate that its reflectivity at normal incidence ranges from 76-80% at 200 nm 
to approximately 90% at 600 nm. 

For tracks in the central region, at polar angles of 65° to 115°, the simple cylindrical 
mirror does not transport sufficient light. Monte Carlo simulations show that the efficiency 
for transporting photons to either end of the tube drops from 30% at 65° to zero at 90°. 
To improve this situation we have developed a specialized reflector for the central region. It 
is designed to redirect a large fraction of the photons such that they exit the region before 
encountering another mirror surface and are then able to reach the end of the module with 
just one or two additional reflections. This reflector has a faceted surface with a sawtooth- 
shaped cross section (Figure ^). Each facet is a planar strip with a width of order 1 mm [Q, 
oriented such that an incoming photon produced by a particle from the interaction point 
is reflected toward the end of the tube, as shown in Figure |[ This faceted mirror system 
extends parallel to the beam line, 0.61 m to either side of the interaction point {z = —0.61 m 
to z = 0.61 m), is about 5 cm wide, and faces the beam line. 

The optimal combination of facet orientations varies with the polar angle of the incident 
particle. Our design incorporates ten different mirror sections, each with area ~ 5 x 5 cm^ 
and a different combination of angles (Table |). The facets reflecting photons toward the 
nearer phototube make an angle a with respect to the z axis, while those reflecting photons 
toward the more distant PMT make an angle f3, as shown in Figure |^. The optimal values 
of a and /? as a function of z depend on the tube length, its radius, the perpendicular 
distance to the interaction point, and the reflectivity as a function of wavelength. A Monte 
Carlo simulation was used to flnd a set of discrete angles for the ten sections spanning 
0.0 m < |;z| < 0.61 m. 

Once a photon reaches the end of the tube, it may need to be funneled to a smaller area 
in order to be delivered to a photocathode; affordable PMT's have typical photocathode 
diameters of 3.8 or 6.4 cm. In addition, a PMT cannot survive a high pressure environment, 
and must therefore be located outside the pressure vessel. The light is delivered to the PMT 
through a quartz window which serves as pressure barrier, UV-transparent window, and light 
funnel. For this test the window was made of Corning excimer grade quartz, in the shape of 
a simple truncated cone, with an inner face slightly smaller than the diameter of the tube 
(9.7 cm), an outer face matching the photocathode (6.4 cm diameter), and a thickness of 
3.8 cm. The conical side wall was coated with aluminum. This window works in two ways: 
because of the change in index of refraction, rays that come in at large angles of incidence 
refract to smaller angles, thereby aligning their trajectories more parallel to the axis of the 
tube; secondly, many of the photons striking the mirrored conical wall are directed into the 
PMT. Monte Carlo studies of photon transport predict that such a window is about 65% 
efficient in transmitting photons from the inner face of the window to the photocathode. 

To minimize transmission loss, the windows of the pressure module and the PMT must 
be coupled to each other by an index matching fluid. Using the monochromator appara- 
tus described in section [11 A| , we compared transmission of light through the window with 
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and without each of four recommended couphng greases, General Electric Viscasil, Oken 
626A, Dow Corning Q2-3067, and Shinetsu. Of the four, only the General Electric Viscasil 
performed well. It improved transmission through the window down to 220 nm and gave 
diminishing efficiency at lower wavelengths, dropping to 80% at 200 nm. The other three 
did not transmit light efficiently below 270nm (less than 40% transmission). 

D. Photomultipier 

In many experimental configurations, the PID detector is located in a strong magnetic 
field. In such cases, conventional photomultiplier tubes cannot function, and it is necessary 
to use a specialized photodetector which is able to operate in such an environment. We 
tested one such PMT, a Hamamatsu R3386 "fine mesh" tube 0, to evaluate its suitability 
for this application. 

For our Cerenkov detector the primary concern is resolution, particularly at low light 
levels; loss of gain may in principle be recovered with additional amplifiers. Signal resolu- 
tion in conventional PMT's is determined mainly by photoelectron counting statistics, as 
well as the quantum efficiency of the photocathode. In fine mesh tubes the resolution is 
complicated because the dynodes are partially transparent (~50%), yielding anomalously 
shaped distributions at low intensities. In fact, a single photoelectron signal will not show 
a well-defined peak, due to this partial transparency |1T0| . 

To evaluate the Hamamatsu R3386, we set up a test system to deliver calibrated light 
pulses to the PMT and then to measure the gain and resolution of the output. A BNC 
light pulser was used to deliver up to several hundred photons, and pulse height spectra 
were recorded using a LeCroy qVt and a Canberra multichannel analyzer. The light output 
of the pulser was calibrated using a Hamamatsu R2059 "photon counting" photomultiplier 
tube. The fractional width of the R2059 was found to be (1.01 ± 0.02)A^~2 , where is the 
number of photoelectrons, in agreement with expectations for a "photon counting" PMT. 

The gain and resolution of the R3386 were then measured in and out of the CLEO II 
solenoid's 1.5 Tesla magnetic field. For this test the axis of the PMT was oriented parallel 
to the field. The high voltage on the R3386 was varied in the range 1200-2000 V at T, 
and 2500-2700 V at 1.5 T, such that all measurements could be carried out within the same 
range of signal heights. The gains and widths at T were extrapolated to the higher voltage 
region for comparison with the data at 1.5 T. Separate measurements in T at the higher 
voltages were used to confirm that the extrapolations were reasonable. For signals of three or 
more photoelectrons delivered to the R3386 in the absence of a magnetic field, a clear peak 
was visible. The fractional width of the signal was found to be equal to (1.10 ± 0.05)A^~2. 
In the 1.5 T field the gain was found to be reduced by a factor of 381 ±65. This corresponds 
to a mean loss of ~31% in multiplication at each of the sixteen stages of the R3386. The 
width of the signal was measured to be (1.25 ± 0.05)A^~2. The responses to a mean signal 
of 10 photoelectrons are shown in Figure H for magnetic fields of T and 1.5 T. 

We conclude that the R3386 is suitable for use in our detector but that it is important 
to take into account the broadening of the signal distribution in predicting identification 
efficiencies. In our detector simulation the photomultiplier was modeled by broadening the 
delivered photon signal by 35%. This is somewhat wider than what we have observed, and 
should therefore give a conservative assessment of feasibility for this detector. 
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E. Electronics 



The electronics to read out this detector needs to satisfy only modest demands. The 
number of channels is low, 200 or less, so that cost is not an issue. Resolutions of ~ 1 ns 
on timing and 5% on pulse height are adequate, and the required dynamic range is three 
decades. 

If the system is to operate in a high magnetic field, preamplification of the phototube 
signal is required, as the gain of the fine mesh PMT may be as low as 10^. For a gain of 10^, 
a single photoelectron will result in signal charge of 1.6 fC on the anode. If this charge is 
spread over 5 ns, it gives a current amplitude of 0.32/iA. A preamplifier gain of 30 mV//xA 
is needed to translate this into a 10 mV signal. The noise at the output should be less than 
1 mV. These criteria are nearly satisfied by the LeCroy TRA4402, which provides a gain 
of 25 mV//iA with an input noise < 65 nA RMS (1.6 mV at the output). There are other 
amplifiers which have input noise current more than a factor of 3 lower than this |]TT| , which 
should easily meet the requirements. 



III. PROTOTYPE TESTS 
A. Optics 

To explore the design parameters and to verify the simulation of this complicated optical 
system, we set up a model of a tubular module with a laser configured to simulate Cerenkov 
radiation. The model tube was half length, 10 cm in diameter, and lined with aluminized 
Mylar D-1 overcoated with MgF2 0. The simulated radiation, an annular "cone" of laser 
light, entered the tube through one of several ports along its length. The light was produced 
with an opening angle of 7.9°, to simulate the radiation from a 1 GeV/c pion. 

To produce the conical distribution, the monochromatic beam from a green HeNe laser 
was first transported through a beam expander with adjustable mask and iris diaphragm 
to generate an annular beam of parallel rays. This beam was then directed onto an axicon, 
which transformed the beam into an annular cone of light. The location for the origin of the 
simulated radiation was adjusted by varying the mask position in the beam expander, and 
the path length over which the radiation was emitted was set by the iris diaphragm. This 
radiation pattern differed from real Cerenkov light in two ways: it did not have the correct 
spectral distribution or polarization. A schematic of the beam and apparatus is shown in 
Figure ^. 

One end of the tube was equipped with an array of photodiodes attached to a movable 
stage that could scan the entire area of the end of the tube. The source intensity was 
measured by shining the cone directly onto the photodiode array. The transport efficiency 
is then the ratio of light intensity detected at the end of the tube to the intensity of the 
source. To correct for temporal fluctuations in laser power, a fraction of the beam was 
extracted via a pellicle beam splitter and monitored during all measurements. To measure 
its intensity in a manner insensitive to spatial fluctuations, the extracted beam was directed 
into an integrating sphere, where it was diffusely scattered from an inner wall coated with 
highly reflective white paint and detected in a photodiode placed out of the beam's path. 
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In our first test, we injected the simulated Cerenkov light at several angles, each at the 
position along the tube corresponding to that expected from particles in an actual colliding 
beam experiment. A 3.0 cm thick plain quartz disk was placed before the detector array to 
simulate a thick transparent endplate. The light intensity was measured over the entire cross 
section of the tube, and the results were compared to expectations from the simulation. The 
measured transport efficiency as a function of angle is shown in Figure |^ with predictions 
calculated by the same Monte Carlo simulation used to design the faceted mirror. The 
agreement between the measured and simulated efficiencies is excellent. It varies from zero 
for a polar angle of 90° to 88% for incidence at 43°, positions corresponding to the center 
and end of the tube, respectively. 

In a detector module the efficiency is further reduced by the fact that even the larger 
PMT photocathode (6.4 cm diameter) covers only ~40% of the total cross sectional area 
of the tube. Thus, with a plain window and the large photocathode, only about 40% of 
the photons reaching the window would strike the photocathode. This yield can be raised 
to around 65% by using the conical window discussed earlier. Data taken with the conical 
window, shown as triangles in Figure ^, are in agreement with this prediction. 

Finally, several 5x5 cm^ prototype sections of the faceted mirror were produced commer- 
cially [Q. A single section, designed for placement at polar angles near 90°, was inserted 
near the far end of the tube. A particle at 90° was simulated and the delivered light measured 
with each of the two windows in place. The transport efficiency was found to be (22±2)% 
for the plain window. The simulation predicts a value of (23±1)%, in good agreement with 
the measurement. These values imply an efficiency of around 8% for delivery to a 6.4 cm 
diameter photocathode. For the conical window the corresponding efficiency was measured 
to be (13±2)%. In a detector module, photons produced by tracks incident at a 90° polar 
angle are reflected equally by the faceted mirrors to both ends of the tube, so the total 
transport efficiency with the conical window will be 26%. 

B. Detector module 

To verify several critical aspects of the design, we constructed one half-length module to 
serve as a prototype (Figure |). Aim long, 10 cm diameter stainless steel tube was lined 
with aluminized Mylar D-1 overcoated with MgF2 0. Two prototype sections of the faceted 
mirror were placed inside the tube near one end, corresponding to the center of a full-length 
module, and the endplate was blackened to eliminate extraneous reflections. Placed at the 
other end of the tube was a Hamamatsu R4143Q photomultiplier with quartz window,^ 
viewing the pressurized gas through the mirrored conical quartz window, described earlier. 
The PMT was optically coupled to the quartz window using General Electric Viscasil. 

To study the response of this prototype module to energetic particles, we constructed a 
cosmic ray telescope, using plastic scintillators to trigger on muons with momentum above 
1 GeV/c. Three small scintillators above the module insured that the cosmic rays passed 



This PMT is not what would be used in a final design, since it does not have the fine-mesh 
dynode structure. 
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directly through the desired region of the gas. An 82 cm thick absorber of lead was placed 
directly below the prototype, and two larger scintillators were placed under the lead. To 
trigger all five scintillators in coincidence, a cosmic ray muon would need not only to pass 
through the prototype at the desired position and angle, but also to have sufficient momen- 
tum to penetrate the lead, at least 1.0 GeV/c. The prototype was supported on a U-channel 
which pivoted about a steel rod such that its angle with respect to the vertical axis and its 
position along the tube axis could be easily modified. A schematic view of the cosmic ray 
test configuration is shown in Figure |[ The photomultiplier signal from the prototype was 
amplified by a Hewlett Packard 5582A linear amplifier and delivered to a Canberra 750Plus 
series pulse height analyzer. The gate for the signal was 4 fis wide, generated in an EG&G 
416A gate generator by a five-fold coincidence of the scintillators. The measured signal was 
calibrated by operating the photomultiplier with its face covered, such that the observed 
signal corresponded to single photoelectron noise. 

The Cerenkov light yields were studied for cosmic rays at many positions along the 
length of the module. In one set of readings, the incidence angle and position were varied 
simultaneously to reproduce those expected of particles entering the actual detector. This 
set consisted of the following combinations of polar angle {6) and distance from the window 
(fe): (a) e = 51.0°, 6z=35 cm, (b) 6 = 58.5°, 5z=50 cm, (c) 6 = 66.8°, 6z=Q5 cm, and (d) 
6 = 73.0°, 6z=75 cm. In the final detector, there would be a faceted mirror module at each 
of the settings (b), (c) and (d), but in this test none was employed. 

These test configurations were duplicated in a Monte Carlo simulation, where particles 
penetrate the prototype and generate Cerenkov photons at the expected rate, distributed at 
the expected angles and energies. Each photon then propagates through the gas volume and 
is absorbed or reflected at each mirrored surface encountered with a wavelength-dependent 
probability based on the measured reflectivity. A photon reaching the front face of the end 
window may be either reflected from the quartz or transmitted. After entering the quartz, a 
photon is reflected back out the front face, absorbed by the mirrored wall, or transmitted to 
the photocathode. Once the photon enters the photocathode, the manufacturer's tabulation 
of its quantum efficiency is used to determine whether a photoelectron is produced. 

The data were compared with the Monte Carlo predictions, and simulation parameters 
were adjusted until agreement was achieved. We found that excellent agreement was ob- 
tained if the reflectivity of the liner in the simulation was adjusted upward by approximately 
1% over all wavelengths. Since the reflectivity used originally as input to the simulation was 
measured on an earlier sample than the liner in the prototype, it is possible that an im- 
provement was made. Figure ^ shows one set of data with its prediction from the tuned 
simulation. 

After verifying agreement between the simulation and data without the faceted mirror, 
we aligned the prototype module horizontally, in an orientation corresponding to a polar 
angle of 90°. For each of several positions, a faceted mirror section was placed such that 
Cerenkov photons would strike it and be directed toward the PMT. The mirror and trigger 
were set up such that cosmic rays entered at average distances of 45 cm, 57 cm, 73 cm, and 
90 cm from the PMT. It is important to note that only the 90 cm position corresponds to 
a physical situation if the module were installed around an interaction region. 

In this configuration, because the area of a faceted mirror section is comparable to that of 
the trigger counters, the raw data include events where most of the light misses the faceted 
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section. To account for such events, each test was repeated with the mirror section covered 
with black paper, and the resulting spectrum was normalized to the number of triggers and 
subtracted directly from the raw data collected with the mirror section uncovered. The 
results are in good agreement with our simulation, if the reflectivity of the mirror sections 
is taken to be about 0.65 times that of the reflective liner. The sections used here were the 
first replicas to be produced from a master. It is expected that with a better understanding 
of the replication process, the reflectivity of the faceted mirrors would increase substantially. 



Results for the 90 cm position are shown in Figure 10. The agreement between data and 



Monte Carlo is excellent, giving us confidence that the simulation is a reasonable representa- 
tion of the prototype and that the critical parameters of this system are well understood. In 
the final system, the light collected from charged particles penetrating near the center of the 
detector module would be approximately doubled, as the module would be approximately 
twice as long and would be instrumented with PMT's at both ends. 



IV. SIMULATION STUDIES FOR CLEO III 

To evaluate the effectiveness of the full detector system for physics results, the Monte 
Carlo simulation described in the previous section was incorporated into the CLEO III 
GEANT-based |]13[ package. It was then used to study two specific physics problems in 



e~^e~ annihilation events in the T region: (1) separation of modes B ^ tttt and B — > Ktt, 
and (2) extraction of — »• ttIu from a large background of — > Klu. 

This simulation includes track curvature due to the magnetic field, interactions at the 
various detector interfaces, splashback from interactions in the Csl electromagnetic calorime- 
ter, and the performance of the fine mesh PMT in a magnetic field. The signal for each 
module is derived by summing the outputs from both PMT's. Since a track almost always 
traverses two modules, the net signal is the sum of outputs from two to four PMT's. The 
noise rate of the fine mesh PMT is of the order of 1 KHz above the 0.5 p.e. level, and the 
event gate in CLEO is estimated to be in the 100 ns range; this implies that there is no 
significant noise contribution when the outputs of more than one PMT are summed. 

First, single vr and K tracks were generated with momenta in the range 0.8-3.0 GeV/c 
and polar angles between 45° and 135° (this angular region is set by the active volume of 
the detector). As an example, we show in Figure |Tl] the expected photoelectron yield for 
pions and kaons of momenta 1.0 and 2.5 GeV/c and polar angle between 85° and 95°. At 
2.5 GeV/c the kaons are slightly above threshold. Beyond the expected signal, it is seen 
that a small fraction of kaons produce a substantial false signal due to an interaction in 
the material before the detector, which may produce a pion or delta ray which is above 
threshold. To be identified as pions (kaons), we require the photoelectron yield to be above 
(below) a certain value. The cut values are different for pions and kaons and are functions 
of the track momentum and polar angle. To reduce the rate of false identifications due to 
interactions, we require alignment between the projected track and the shower centroid as 
measured in the calorimeter. The degree of alignment required is in the range 1.7°-2.3° 
and depends primarily on the track momentum. It is chosen such that 80% of pions pass 
it. In Table |I| we show the net pion efficiency, e,r, and kaon fake rate, fx, where the cuts 
have been optimized by maximizing eT^/fx, which is proportional to the ratio of signal to 
background (S/B). At momenta 2.5-2.8 GeV/c, S/B ranges from about 30 to 60, depending 
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on the angle of the track; this translates roughly to 2.5a tt/K separation. Combining this 
information with measurements of specific ionization in the drift chamber will increase the 
7t/K separation to above 3a . 



Using the criteria developed from simulations of single vr and K tracks, we studied some 
specific decay modes in fully simulated events. The efficiencies presented below are due only 
to particle identification, since all other efficiencies, such as geometric cuts, kinematic cuts, 
lepton identification, and isolation of tracks to eliminate multiple hits, cancel to first order 
for this evaluation. 

We have studied the two decays B ^ and B — > Kn, where each mode poses the 
most serious background for the other and the two branching fractions are thought to be 
comparable. If we choose pulse height cuts by maximizing S/B, we find that the PID 
efficiency for detecting the Kn final state is (23.2 ± 1.0)%, whereas the background from 
the TTTT state is (0.3 ± 0.1)%. If instead, we wish to detect B vrvr, then the efficiency is 
(22.7 ± 1.0)% and the background from the Kit final state is (1.2 ±0.1)%. These efficiencies 
imply a rejection ratio of at least 1/20, which corresponds to a separation of 2.5a. The S/B 
for the two modes are different because kaons have a long tail at high photoelectron yields, 
while pions do not have a similar tail at low yields. Further separation between the two 
samples can be achieved by using dE/dx information |T^ from the drift chamber and the 
kinematic quantity AE = Eheam — Eb- 

We also studied the decay n^l^u, a Cabbibo-suppressed decay where there is 

a large background from the unsuppressed decay D° — K^l^u. In this case, cuts were 
optimized by maximizing S^/ (S+B), where the vr mode was assumed to have a rate 15 times 
smaller than the K mode. We selected D mesons produced with momenta above 1.9 GeV/c, 
and with n or K polar angle between 45° and 135°. The resulting efficiencies for detecting the 
pion in ti^I^u and fake rates from kaons in D° K~l^v are presented in Table |T|. 

These results indicate a S/B of about 20. It is clear that, even after allowing for the higher 
rate for —>■ K~l^v, the final sample is enriched with the signal decay. This is much 



better than the current CLEO II performance 15 , where the final sample is dominated by 



background from K l^p events. Additional signal enhancement may be achieved through 



kinematic requirements and measurements of specific ionization WA 



V. CONCLUSION 

It has been demonstrated that a simple, compact gas threshold Cerenkov counter which 
utilizes commercially produced photodetectors can provide hadron identification in the mo- 
mentum region 1-3 GeV/c. The design described here can be implemented as is, with the 
development of a mechanical design. Possible improvements requiring some further devel- 
opment include better quality refiective liners and faceted mirrors, optimized design of the 
quartz window, and use of low-Z materials such as carbon fiber or beryllium in the tubes to 
reduce the thickness of material. In addition, continuing improvements in photodetection 
technology are beginning to yield devices with better operating characteristics, e.g., higher 
quantum efficiency, than the fine mesh PMT studied here ||16|| . 

Although this detector was developed for implementation in the CLEO III detector, it 
could be used in a much broader set of applications. Because of its speed and the low 
number of readout channels, it can be used in most settings with moderately high rate and 
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low particle multiplicities. By using multiple layers and varying gas pressures, this design 
can be applied to all particles with momenta above 1 GeV/c, up to 10 GeV and above. These 
are difficult regions in which to achieve hadron identification, and the detector presented 
here does so simply and at moderate cost. 
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FIG. 1. Transmission curves for two grades of SFq at a pressure of 20 atm. Because of low 
light yields, each point has an error of ~4%, which is shown for only two data points. 
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FIG. 2. Configuration of detector modules in a coUiding beam experiment. 
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FIG. 3. Cross section of faceted mirror. The nearer photomultiplier is toward the right in the 
figure. 
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FIG. 4. Process for transport of photon from center of detector module to photomultiplier. 
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FIG. 5. Signal distributions for R3386 fine mesh phototube in T and 1.5 T magnetic fields, 
for light pulses producing an average of 10 photoeiectrons. Both signals are normalized to a mean 
of 10 photoeiectrons. 
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FIG. 6. Schematic diagram of optical test configuration. 
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FIG. 7. Photon transport efficiencies for green light with plain cylindrical reflector at large track 
polar angles. Results for light collected over the entire area of a quartz disk end window (10 cm 
diameter) show good agreement between data (squares) and Monte Carlo prediction (circles). In 
a real detector the photosensitive area will be no larger than 6.4 cm diameter, so the yields for the 
plain disk will be 40% of those shown. Also shown is the result from data taken using a mirrored 
conical quartz window (triangles), where light is collected over a 6.4 cm diameter area. 
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FIG. 9. Photoelectron yields for cosmic rays entering the prototype module at angles which 
correspond to a polar angle of (a) 51°, and (b) 58.5°. The points with error bars are data, and the 
solid histogram is the Monte Carlo prediction. In the design for CLEO III, the light from a track 
entering at an angle of 58° would encounter a faceted mirror. In this test there was no faceted 
mirror, thus giving a lower light yield. 
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FIG. 10. Photoelectron yield in the prototype. Cosmic rays enter at an angle which corresponds 
to a polar angle of 90°, at a distance of 90 cm from the photo multiplier tube, and strike a 
faceted mirror. The points with error bars are data, and the solid histogram is the Monte Carlo 
prediction. In the final system, this yield would be approximately doubled, as the module would 
be approximately twice as long and instrumented with PMT's at both ends. 
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FIG. 11. Monte Carlo prediction of signal distributions for pions and kaons with momenta of 
1.0 GeV/c and 2.5 GeV/c incident on the detector with polar angle between 85° and 95°. 



22 



TABLES 



TABLE I. The orientations of the mirror facets, a 


and [3 are defined in 


Figure y 


\z\ interval (m) polar 


angle (degrees) 


a (degrees) 


j3 (degrees) 


0.00 - 0.05 


88.4 




37.5 


37.5 


0.05 - 0.11 


84.9 




32.0 


41.0 


0.11 - 0.15 


81.8 




29.0 


42.0 


0.15 - 0.21 


78.7 




26. U 


AO r\ 
43. U 


0.21 - 0.26 


75.4 




23.0 


44.0 


0.26 - 0.32 


72.1 




20.0 


45.0 


0.32 - 0.38 


68.7 




18.0 


46.0 


0.38 - 0.45 


65.2 




16.0 


47.0 


0.45 - 0.53 


61.4 




14.0 


48 n 


0.53 - 0.61 


57.7 




12.0 


49.0 


TABLE IL Pion efficiencies e 


and kaon fake rates fx 


in percent, for sing 




different polar angles. Cuts have been optimized by maximizing S/B. 




Momentum (GeV/c) 


50 ±5° 




70 ±5° 


90 ±5° 




Ctt Uk) 




evr Uk) 


evr Uk) 


0.8 


54.0 (0.65) 




70.8 (0.90) 


62.8 (1.35) 


1.0 


54.0 (0.25) 




72.1 (1.40) 


61. y (l.lUj 


1.2 


59.6 (0.20) 




71.2 (0.80) 


62.4 (0.80) 


1.6 


59.9 (0.87) 




70.6 (1.15) 


62.7 (1.20) 


2.0 


54.6 (0.65) 




70.3 (1.50) 


Di.i (U.OOj 


2.5 


37.6 (1.05) 




71.6 (1.30) 


61.2 (0.95) 


2.8 


36.3 (1.10) 




42.2 (1.55) 


28.3 (0.90) 


3.0 


36.5 (2.15) 




32.1 (1.30) 


33.5 (3.50) 


TABLE III. Efficiencies for 


TT l^u, fake rates from K 


t v 


Track momentum (GeV/c) 




Efficiency, nlv 


Fake rate, Klv 


0.8-1.0 




70.7% 




4 6% 


1.0-1.2 




67.3% 




4.5% 


1.2-1.6 




71.3% 




2.6% 


1.6-2.0 




74.3% 




4.1% 


2.0-2.8 




71.1% 




2.9% 


2.8-3.0 




46.8% 




4.4% 


0.8-3.0 




(70.5 ± 1.3)% 


(3.6 ± 0.3)% 



23 



